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A supersymmetric composite model of color superconductivity is proposed. Quarks and diquarks 
are dynamically generated as composite fields by a newly introduced strong gauge dynamics. It is 
shown that the condensation of the scalar component of the diquark supermultiplet occurs when the 
chemical potential becomes larger than some critical value. We believe that the model well captures 
aspects of the diquark condensate behavior and helps our understanding of the diquark dynamics in 
real QCD. The results obtained here might be useful when we consider a theory composed of quarks 
and diquarks. 
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In the past few years, interesting properties of quark 
matter with nonzero temperature and baryon density 
such as the phase diagram in Fig. ^ have been exten- 
sively studied. In particular, color superconductivity Q 
is one of the most promising possiblities to occur in such 
a system at low temperature because there exists an at- 
tractive interaction between quarks on the Fermi surface 
through the antisymmetric color 3 gluon exchange and 
it results in the quark-quark (not quark-antiquark) pair- 
ing, so-called diquark condensate. So far, a lot of people 
have tried to investigate this phenomena using different 
methods. At asymptotically high density, one can per- 
form the weak coupling analysis of QCD and compute 
the superconducting gap [2| • On the other hand, at lower 
densities where the quark chemical potential // is of or- 
der 500 MeV, since QCD itself is not tractable in this 
strong coupling regime and in addition there is a norto- 
rious sign problem on the lattice Monte Carlo simulation 
with nonzero [A, only some model studies (for instance, 
the Nambu- Jona-Lasinio model) which mimic crucial fea- 
tures of QCD such as chiral symmetry have been done so 
far [3 . According to these model studies, a phase transi- 
tion between hadronic phase and color superconducting 
one is suggested to happen. Combining this with the re- 
sult of the weak coupling analysis, we expect a similar 
phase structure is realized in dense QCD. 

In the light of these current situations, the purpose of 
this paper is to propose some new way of understanding 
a system of dense quark matter based on supersymmetric 
(SUSY) QCD. In SUSY gauge theories, some exact non- 
perturbative results, which are powerful tools to study 
the strong coupling dynamics, have been already known 
4] . It might be therefore interesting to investigate color 
superconductivity using SUSY gauge dynamics and con- 
sider whether we can obtain some insight into real QCD. 



In symmetry breaking pattern was studied by us- 
ing the exact results of SUSY QCD with nonzero chem- 
ical potetnial and was compared to that obtained from 
the analysis via nonsupersymmetric QCD One of 
the important observations is that the chemical poten- 



FIG. 1: (Conjectured) phase diagram of hot and dense quark 
matter. 



tial can be incorpolated as the time component of a ficti- 
tious gauge field of the U(1)b baryon number symmetry 
at zero temperature, which leads to a tachyonic SUSY 
breaking scalar mass. However, gauge variant quantities 
such as the diquark degrees of freedom was not treated 
in their model. In order to see what happens to the sys- 
tem at large chemical potential, however, as has been 
already mentioned before, it is indispensable to include 
that degrees of freedom. Therefore in this paper we try 
to extend the work of Q so as to involve gauge variant 
quantities at an intermediate energy range. 

Along these line of thought, we propose a supersym- 
metric composite model of color superconductivity, in 
which quarks and diquarks appear as massless compos- 
ites at low energy by a newly introduced strong cou- 
pling gauge dynamics, not by QCD dynamics. We find 
a certain parameter region where the scalar component 
of diquark supermultiplet condensation occurs when the 
chemical potetnial gets larger than some critical value. 
Although our model is not fully realistic in the sense 
that not diquarks themselves but the scalar component 
of the diquarks supermultiplet condense, nevertheless we 
believe that our model well captures some important as- 
pects of the diquark condensate and helps our under- 
standing for the color superconductivity in real QCD. 
The results obtained here may give an interesting insight 



on the phase structure at the intermediate region of the 
quark chemical potential. 

Let us explain our model, which is based on a M = 1 
SUSY SO(N = N f + 4) gauge theory with N f vector 
representations Q. Its non-Abelian global symmetry 
SU(N f ) is extended to SU(3) C x SU(N f ) L x SU(N f ) R 
where SU(3)c (usual color symmetry) gauge theory is as- 
sumed to be weakly gauged compared with SO(N) gauge 
theory, l.e.Aso(N) > ^SU(3) C f° r the dynamical scales of 
each gauge group. Matter content is summarized below. 



Q = (□,□,□, l)i, 1,1, 
Q = (□,□,].,□ )-i,i,i, 

X = (□ , 1, l,l)o,-6JV/,3— JV 
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where the representations in the parenthesis are transfor- 
mation properties under the group SO(N) x SU(3)c X 
SU(N f ) L x SU(N f ) R , where N = 6N f + 5 in the present 
case. The numbers in the subscripts are charges for 
nonanomalous U(l) global symmetries U(1)b X U(1)a X 
U(1)r, which each U(l) symmetries are linear combina- 
tions of the original anomalous U(l) symmetries. U(1)b 
is a baryon number symmetry which plays an impotant 
role for considering the chemical potential effects. U(1)a 
is a non-R symmetry and U(1) R is an R-symmetry. 

SO(N) gauge theory under consideration is asymptot- 
ically free and known to be in a confining phase at the in- 
frared (IR) 0). At the scale A. SO (n), the theory becomes 
strongly coupled and SO(N) gauge invariant composite 
fields appear as massless degrees of freedom in the low 
energy effective theory. 
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where the representations in the parenthesis are those 
under the group SU(3) C x SU(N f ) L x SU(N f ) R . The 
numbers in the subscripts are charges for nonanoma- 
lous U(l) symmetries U(1) B x U(1) A x U(l) R . Note 
that Q 2 (Q 2 ) has symmetric (its conjugate) and 
anti-symmetric (its conjugate) representations under 
SU(3) C x SU(N f ) L x SU{N f ) R because SO(N) in- 
dices are contracted symmetrically and the superfields 
are bosonic. The anti-symmetric ones correspond to 
"diquark" superfield responsible for the condensation 
when the chemical potential becomes larger than some 
critical value. QX(QX) correspond to usual quarks 
(anti-quarks) superfields. Thus, quarks and diquarks 
coexist as massless composites in the low energy. As 
a nontrivial check that composite fields are 
appropriate massless degrees of freedom, we can easily 



show that the 't Hooft anomaly matching conditions 
for [SU(N f ) LtR ] 3 , [SU(N f ) L . R } 2 U(l) B .A.R, U(l) B , A ,R, 
[U(1)b,a,r] 3 , [U(1)b,a] 2 U(1) r , [U(1)am] 2 U(1) b , 
[U(1)b,r} 2 U(1)a at the origin of the moduli space 
are satisfied between elementary fields Q,Q,X and all 
composite fields (Q}-©. 

The low energy effective superpotential is generated by 
the gaugino condensation in the unbroken gauge group 
SO(N) - 50(4) ~ SU(2) L x SU(2) R ; 



W cB = 2{e L + e R ) 



( A 6A'/+4 \ 
/ 1V SO(N) \ 

\[det(QQ)]X j 



(10) 



where £l,r. — il are phase factors reflecting the number 
of SUSY vacua suggested from Witten index || . 

Since our interest is whether the condensation of the 
scalar component of the diquark supermultiplet occurs 
or not as the chemical potential changes, we need to es- 
timate the soft scalar mass squareds for composite fields, 
whose sign indicate whether composite fields develop vac- 
uum expectation values (VEVs) or not. In softly bro- 
ken SUSY gauge theory, it is well known that soft scalar 
masses in the IR region can be derived from those in the 
ultraviolet (UV) region by the procedure in Ref. @. The 
effective Kahler potential for composite fields is fixed by 
symmetries and the renormalization group (RG) invari- 
ance, 



C (Q 2 ) j w ) e (Q ) 



+ c (02) ^p(Q 2 )t e -^(Q 2 ) 

+ C( * 2 )fe(Y 2 )t(Y 2 ) 

+ c (QX) ^f^(QX)^e VB (QX) 
Z, 

+ C (QX) j 
Z, 



(11) 



where overall coefficients c's are of order 0(1) unknown 
constants. The exponential factors for QCD are sup- 
pressed. Vb is a background vector superfield U{1) B 
with a VEV (V B ) = fo»B(Ap), (A^) = (ji/g B , 0,0,0). 
cjb is a gauge coupling constant and \i is a chemical po- 
tential, which breaks SUSY explicitly. Therefore, we 
assume /1 <C Aso(N) so that we can make use of ex- 
act results of the SUSY gauge theory. This VEV of 
the background vector field provides additional tachyonic 
soft SUSY breaking scalar mass squareds as discussed in 
|5|- Wave function renomalization constants Zi are pro- 
moted to a superfield Zj 



Z, = Zi [1 - 6 2 9 2 m 2 soit ] 



(12) 
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where m so ft is a soft SUSY breaking scalar mass in the 
UV and taken to be universal. The quantity / is a spu- 
rious U(l) symmetry and the RG invariant superfield, 



I = k\z 2T ' b °k h 



(13) 



where T is the total Dynkin index of the mat- 
ter fields, 60 is the 1-loop beta function coeffi- 
cient and A h = fi uv exp[-&TT 2 S(fj, uv )/b ], S(fi uv ) = 
p- (l + # 2 ^) (m\ : gaugino mass). Note that a spu- 
rious U(l) transformations are given by 



e A+Ai Z ri 
T 



SVuv) -> S(n vv ) - —A 



(14) 
(15) 
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FIG. 2: Relation among various scales are displayed. Hor- 
izontal axis means the energy scale. The possible range of 
the chemical potential is below v/ m 2 + . There is no 
condensation when the chemical potential is between QCD 
scale Agu(3) and the critical chemical potential n*. The con- 
densation of the scalar component of the diquark supermul- 
ti plet occurs when the chemical potential is between fi, and 
sj m 2 + m 2 o{t . UV theory becomes unstable if the chemical 
potetntial is beyond the scale W rn 2 + m 2 oft . 



where A is a chiral superfield. 

Now, the soft SUSY breaking scalar masses for com- 
posites with the canonical kinetic term are obtained by 
taking 6 2 9 2 terms in Eg. (|H)> 



term l)19|l can be rewritten in the IR as the linear term 
of composite, 



TO (Q 2 ) - m (Q 2 ) 
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(16) 



W = m(QQ), 



(20) 



"4ft > 0, (17) 
TO 2 oft - 7 /A (18) 



For the case cl = —£r with vanishing superpotetnial in 
l|l()|l. the scalar potetntial consists of only the soft scalar 
masses Q16 [I -I|1H |1 . We then immediately find from i|16|) 

that m 2 Q2) = rh 2 Q 2) < for fi > /Lt*(= \J 2 (3N^+2) m ^ > 
"^soft)- This means that (Q 2 ), (Q 2 ) ^ in that range of 
the chemical potential |Tlj . 

Furthermore, if we take into account the most at- 
tractive channel hypothesis fl^ ]. anti-symmetric part of 

Q 2 (H ,0 , 1)2,2,2, Q 2 (n , 1, )-2,2,2 are likely to have 
VEVs since the force acting on Q 2 or Q 2 by one SO(N) 
gauge boson exchange is attractive. On the other hand, 
the force is replusive in the symmetric case. This is our 
main result that we wish to show. 

We also note that /i* > m so f t implies /i > m so f t for 
the condensation of the scalar component of the diquark 
supermultiplet to occur. This leads to UV unstable the- 
ory since the soft SUSY breaking scalar mass squared in 
the UV becomes negative in the presence of the chemical 
potential, ttiq = 
mass terms 



■fj? < 0. We therefore add SUSY 



U = in,, (J, (J,. 



(19) 



where the mass is taken to be flavor diagonal to preserve 
the most global symmetry. For the UV theory to be 
stable, y/ m 2 + m 2 o{t > \x is required. The SUSY mass 



which breaks SUSY spontaneously. We also have to im- 
pose the condition m Aso(N) to use the exact results 
of SUSY gauge theory reliably. 

Combining these facts, we can obtain the diquark su- 
permultiplet the scalar component of the diquark super- 
multiplet condenses at some critical value of chemical 
potenial /i*. We cannot determine exactly where the 
scalar component of the diquark supermultiplet stabi- 
lize, but in a region where the field VEV is much larger 
than Aso(N), elementary fields are appropriate variables, 
we know that the potential is stabilized by SUSY mass 
terms. This implies that the VEVs of the scalar compo- 
nent of the diquark supermultiplet should be stabilized 
at certain value. 

One may find from l|18|) that composite squarks de- 
v elop VEVs (QX), (QX) ^ for fi > 2/^ if 2(jl» < 
yj m 2 + T^goft ■ However, the scalar component of the di- 
quark supermultiplet condenses as far as the chemical po- 
tential is larger than the critical chemical potential value 
/i*. In order to analyze the theory in such a case, we 
have to expand our theory around the VEV of the scalar 
component of the diquark supermultiplet condensation. 
Therefore, the above description of the squark condensa- 
tion is left untouched at the present stage. 

The above argument for the behavior of the scalar com- 
ponent of the diquark supermultiplet or squark supermul- 
tiplet condensation is valid for the vanishing superpoten- 
tial with €1, — —6r in i|l(J|) . For the case with = £r, 
on the other hand, it is found that the scalar potenial is 
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very complicated, 
' d 2 K cS 



V = 



dW cS \ ( dW r 
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where $i denote the scalar component of composite su- 
perfields, K e g is given by l(TT|l and W c s is the sum of 
the superpotential (|10|l rewritten in terms of composites 
and (|20|l . Therefore, we give here a qualitative discussion 
on the scalar potential behavior instead of performing an 
explicit minimization of the scalar potential. Note that 
F-term contributions to the scalar potential in the first 
line of 121|) have a runaway behavior, which make the 
fields VEV away from the origin. For \l = 0, all SUSY 
breaking scalar mass squareds are positive, which set the 
fields VEV at the origin. Therefore all composites are ex- 
pected to develop nonvanishing VEVs by balancing terms 
between the runaway potential and the SUSY breaking 
scalar mass terms. Even if we take into account that 
the scalar diquark mass squareds become negative for 
//>//*, qualitative features of phase transition remains 
unchanged. In any case, the case of nonzero superpoten- 
tial with €l = 6r in (|TU1) is irrelevant to the phase of color 
superconductivity of our interest. Even if we compare the 
vacuum energy in both cases, the case with vanishing su- 
perpotential seems to be energetically favored. 

In summary, motivated by the work of we have 
tried to construct a toy model where gauge noninvariant 
operators are taken into account in SUSY gauge theories. 
We have proposed a SUSY composite model of color su- 
perconductivity, which is based on an SO(N) gauge the- 
ory with (N — 4) vector representations pj. Our model 
is in a confining phase for SO(N) gauge dymanics at low 
energies, in which quarks and diquarks are generated dy- 
namically as composite fields satisfying anomaly match- 
ing conditions. We have shown that the scalar component 
of the diquark supermultiplet condensate occurs when 
the chemical potential becomes larger than some criti- 
cal value ju* = y / 2 (3 j jv / + 7 2) m soft- 



Our model is valid in 



the parameter region Agu( 3 ) < \i < yj m 2 + m 2 oft , where 
the upper bound is required for the theory to be stable 
in the UV and the lower bound implies that we con- 
sider the theory where quarks are deconfined. Although 
the model is not fully realistic in that the scalar compo- 
nent of diquark supermultiplet (not diquarks themselves) 
condense, we believe that it well captures some impor- 
tant aspects of the diquark condensation behavior and 
helps our understanding for the color superconductivity 
in real QCD. If there is a certain intermediate region 
of the chemical potential where quarks are deconfined 



but not superconducting yet, owing to the strong quark- 
quark correlation, the system may be well described by 
a compositon of quarks and diquarks. Then the analysis 
performed in this paper will help us with comprehending 
the behavior of such a system. 

As future directions, it is interesting to extend our 
analysis to other flavor cases with various phases other 
than the confining phase. In particular, it might be possi- 
ble to obtain better and more realistic understanding for 
the diquark condensation behavior by exploiting Seiberg 
dual magnetic description In order to fully under- 
stand the phase structure of QCD, it is necessary to take 
into account the finite temperature effects. It is there- 
fore indispensable to consider our model extended to five 
dimensional spacetime compactified on S , and then to 
study the scalar component of the diquark supermulti- 
plet condensation behavior on the temperature-chemical 
potential plane as shown in Fig^ 
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